Lattice dynamics of FeSb 2 
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The lattice dynamics of FeSb2 is investigated by the first-principles DFT calculations and 
Raman spectroscopy. All Raman and infrared active phonon modes are properly assigned. 
The calculated and measured phonon energies are in good agreement. We have observed 
strong mixing of the A g symmetry modes, with the intensity exchange in the temperature 
range between 210 K and 260 K. The A 9 modes repulsion increases by doping FeSb2 with 
Co, with no signatures of the electron-phonon interaction for these modes. 
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I. INTRODUCTION 

FeSb2 is a strongly correlated narrow-gap semiconductor which has recently attracted a lot of 
attention due to its unusual thermoelectric properties. [l-|fi|It was shown that FeSb2 has colossal 
thermopower S at 10 K (range from 1 mV/K to 45 mV/K aQ]) and the largest power factor S 2 a 



ever reported. The phonon contribution to S remains controversial. [10] Also, the thermal 

conductivity k of FeSb2 is relatively high and around 10 K is dominated by the phonons.[5] Con- 
sequently, full knowledge of FeSb2 lattice dynamics is necessary in order to understand the low 
temperature transport and thermodynamic properties of this material. 

The infrared active phonon frequencies of FeSb2 were obtained from the polarized far-infrared 
reflectivity spectra. [ll] From E||b polarized reflectivity measurements on (102) plane of FeSb2 
single crystal, Perucchi et al. observed four modes at 106.4, 231, 257 and 271 cm -1 at 10 K (factor 
group analysis predicts 3-£>2m modes for this polarization configuration). For E_Lb polarization, 
both B\ u and E>3 U symmetry modes can be observed from (102) plane. Three (of four) modes 
at 121, 216 and 261.4 cm -1 are observed for this polarization. Raman scattering measurements 
on FeSb2 were published in Refs. [l2-15. Lutz and MullerfljJ] observed two Raman active modes 



at about 175 and 154 cm 1 on hot-pressed samples, and assigned them as the A„ symmetry 



modes. In contrast, Racu et aZ. [13[| observed three Raman modes at about 150, 157 and 180 
cm -1 using polarized Raman scattering measurements on FeSb2 single crystals and assigned them 
as the Bi 9 , A g and Bi ff symmetry modes, respectively. Finally, all six Raman active modes of 
FeSb2 predicted by the factor group analysis (2A 3 +2Bi 9 +B-2 ff +B3 9 ) were observed in Ref. Il4 . 



Polarized Raman scattering spectra of the Fei_ a; M. r Sb2 (M=Cr,Co) single crystals was studied in 



Ref. [15j. The linewidths and energies of the Raman modes were analyzed as a function of doping 
x and temperature. Strong electron-phonon interaction, observed for the Bi 9 symmetry mode of 
pure FeSb2, produces significant mode asymmetry. The electron-phonon interaction is drastically 
reduced with increasing concentration of Co and Cr in Fei_ a; (Co,Cr) a ;Sb2. The mixing of the A g 
symmetry phonon modes has been observed both in pure and Cr-doped samples. {3] 

In this paper we report ab initio study of the lattice dynamics of FeSb2- The calculated phonon 
energies in the V point are in good agreement with experimental data. Phonon density of state 
show a gap at about 175 cm -1 , which divides a low frequency region where vibration modes are 
mostly Raman active from a high frequency region where only infrared active modes appear. The 
calculated phonon dispersions for two A g symmetry modes indicates strong mode mixing. This 
is indeed observed in our polarized Raman scattering spectra. The A g mode intensity exchange 
in the temperature range between 210 K and 260 K agrees well with theoretical calculations, 
excluding any additional temperature dependent electron phonon coupling for these modes. The 
mode repulsion increases with Co doping. 



II. EXPERIMENT 



Single crystals of FeSb2 and Feo.75Coo.25Sb2 were grown by the self-flux method and character- 
ized as described elsewhere. [l| The Raman scattering measurements were performed using Jobin 
Yvon T64000 Raman system in micro-Raman configuration. The 514.5 nm line of an Ar + /Kr + 
mixed gas laser was used as an excitation source. Focusing of the laser beam was realized with a 
long distance microscope objective (magnification 50x). We have found that laser power level of 
0.02 mW on the sample is sufficient to obtain Raman signal and, except signal to noise ratio, no 
changes of the spectra were observed as a consequence of laser heating by further lowering laser 
power. The corresponding excitation power density was less then 0.1 kW/cm 2 . Low temperature 
measurements were performed between 15 K and 300 K using KONTI CryoVac continuous Helium 
flow cryostat with 0.5 mm thick window. Raman scattering measurements of pure and Co doped 
FeSb2 samples were performed using the (101) oriented FeSb2 samples. Selection rules for parallel 
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FIG. 1. (Color online) Atomic displacement patterns for the vibrational modes 
arrows are proportional to square root of the vibration amplitudes. 



of FeSb 2 . The length of the 
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and crossed polarization from the (101) 
assignment have been presented in Ref. 



alane of the orthorhombic crystal symmetry and the mode 
14. 



III. NUMERICAL METHOD 



FeSb2 crystallizes in the orthorhombic marcasite-type structure of the centrosymetric Pnnm 



(D^) space group, with two formula units (Z=2) per unit cell 



Basic structural unit is 



built up of Fe atoms surrounded by deformed Sb octahedra. These structural units are corner 
sharing in the (ab) plane and edge sharing along the c-axis. Two Fe atoms are in (2a) Wyckoff 
positions at (0,0,0) and four Sb atoms are in (4g) Wyckoff positions at (0,u,v) of the Pnnm space 
group. Our density functional theory (DFT) calculations are performed within generalized gradi- 
ent approximation (GGA) with PW91 exchange-correlation functional which is used to calculate 



ultra-soft pseudopotentials, 



IS] as implemented in the QUANTUM ESPRESSO package. [IJ] Iron 



(antimony) pseudopotential takes into account 3s 2 3p 6 4s 2 3d 6 (4d 10 5s 2 5p 3 ) electron states for the 
valence electrons. The Brillouin zone was sampled with an 8x8x8 Monkhorst-Pack k-space mesh 
and with the Marzari-Vanderbilt cold smearing (0.005Ry).[3] The obtained optimized structural 
parameters are a=5.859 A, 6=6.583 A, c=3.812 A, u=0.1882, and w=0.3554, which are in good 
agreement with the experiment. Our band structure calculations agree well with the previously 



reported. 
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IV. RESULTS AND DISCUSSION 



23] 



The lattice dynamics is investigated by the density functional perturbation theory (DFPT) 
within the theory of linear response. This method includes calculations of charge response to 
the lattice distortions (allowed by the symmetry operations) for the specified vectors in the first 
Brillouin zone. Calculations start from the previously calculated ground state atomic and electronic 
configuration and continue with the self-consistent calculations of the charge response for each 
different displacement. The normal modes of the optical active phonons (in the T point) are given 
in Fig. [TJ Because Fe ions are located in the center of inversion of Pnnm space group, they do not 
contribute to the Raman scattering process, i.e. the Raman modes of FeSb2 originate only from 
the Sb atoms vibrations, in a manner illustrated in Fig. [TJ In the case of infrared active modes, 
both the Fe and Sb atoms contribute to the normal modes, see Fig{TJ 

In order to obtain the phonon dispersion curves, we have calculated the phonon frequencies at 
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TABLE I. Raman and infra-red active mode energies (in cm 1 ) of FeSb2 single crystal. 



Exp.[14j n(o) 



Activity 



Symmetry 






Calculation 




A 1 


150.7 


160.3 


146.8 
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153.6 


164.4 


161.2 


R 




154.3 


164.6 


173.3 


R 




173.9 


190.4 


184.3 
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B 2g 


90.4 




89.4 
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B 3s 


151.7 




95.6 
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Calculation 



Bi„ 



R 3 



B 3u 



R 3 



A 1 

A ol 



195 
106.4 
231.0 
257.0 
271.0 
121.0 
216.0 
261.4 



177.4 
108.1 
210.7 
259.6 

122.8 
219.2 
236.6 
84.9 
195.2 



IR 

IR 

IR 

IR 

IR 

IR 

IR 

IR 
Silent 
Silent 



4x4x4 Monkhorst-Pack (/-points mesh and interpolated along the chosen path. Figure (2ja) shows 
the calculated phonon dispersions, whereas Fig. &b) represents the phonon density of states of 
FeSb2 • It is interesting to note that there is a frequency gap in the phonon dispersion of FeSb2 at 
about 175 cm -1 . The lower frequency range is dominated by Sb-atoms vibrations (mostly Raman 
active vibrations), whereas the Fe atoms vibrate at frequencies higher than 175 cm -1 . These modes 
are only infrared active. 

The phonon density of states peaked structure between 50 an 90 cm -1 correspond to the low 
frequency acoustic modes associated with the low-lying B2 9 Raman active mode, which calculated 
frequency is at 89.4 cm -1 . Sharp peaks in the phonon density of states above 90 cm" 1 come from 
the flat regions of dispersion curves of corresponding Raman (below 175 cm -1 ) and infrared (above 
175 cm -1 ) modes. 

The lattice dynamics calculations allow us to assign the infrared active modes, experimentally 
Lj. The assignment of the infrared active modes is done according to the mode 



observed in Ref. 
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energy and symmetry. As we have already mentioned, for E\\b (E>3 U symmetry modes) four modes 



are observed 



111 ] instead of three. We believe that the appearance of two modes at about 257 and 



271 cm 1 instead of a single frequency mode is the consequence of splitting of relatively broad 
oscillator (which calculated TO frequency is 259.6 cm" 1 ), due to anharmonicity effects. [241] The 
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k-point grid Dos (states/cm 1 unit cell) 

FIG. 2. Phonon dispersion (a) and phonon density of states (b) for FeSb2. Triangles and circles represent 
experimentally observed Raman and infrared mode energies. Inset: Dispersion curves of A g phonon modes 
along the T — Y and T — X directions. 



E>i u infrared active mode of FeSb2 is recently observed at 195 cm 1 in Ref. 25. The frequencies and 
assignment of all infrared active modes are given in Table HI 

The DFT calculations are performed at zero temperature and should be matched with the 
phonon energies at zero temperature. For this purpose, we have analyzed the change of the Raman 
mode energy and linewidth with temperature, induced by anharmonicity effect. 

The influence of the anharmonic effects on the Raman mode energy can be taken into account 
via three- and four-phonon processes by applying the Klemens's ansatz:^. 



fi(T) = O - A^(T) - A< 4 '(T), 

2 



A (3 \T) =c[ 1 + 

a^\t) = d(i + 



1 



(1) 



+ 



e v - 1 {ev - l) 2 



where VLq is the temperature independent contributions to the Raman mode energy, C (D) is the 
three (four)-phonon anharmonic constant, x = fiVL$/2kBT and y = hQ^/SksT. 

There are two main contributions to the phonon linewidth: (i) anharmonic decay of the phonon, 
and (ii) perturbation of the translational symmetry of the crystal by the presence of impurities 
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FIG. 3. (Color online) Ti\ g mode wavenumber as a function of temperature. Solid line represents theoretical 
extrapolation by using Eq. ([1]), whereas dashed line represent theoretical extrapolation obtained by omitting 
four-phonon contribution in Eq. (p}. Inset: Theoretical calculation of FWHM obtained by using Eq. ([2]) (solid 
line) and by omitting four-phonon contribution in Eq. (dashed line). 

and defects. Having this in mind, the phonon linewidth can be described with: 

r(T) = r + r^{T) + r^(T), 
r (3) (T) = ^(i + -^) (2) 

r(4) ^< 1 + ^ + (^) 

where Tq is the temperature independent linewidth, which originates mainly from (ii), A (B) 
is the three (four)-phonon anharmonic constant. Analysis of energy and FWHM (full width at 
half maximum) vs. temperature for the B\ g mode is presented in Fig. [3j Because anharmonicity 
constants ratio B/A and D/C is very small, see Fig. [3], the contribution of the four-phonon processes 
is small compared to that of the three-phonon processes. The obtained value of £1(0) =164.6 
cm -1 for this mode at zero temperature is in good agreement with the DFT results. Similar 
analysis have been performed for symmetry mode giving the value of f2(0) =190.4 cm -1 at 
zero temperature, |15(| which is in rather good agrement with our calculations. 

The calculated energy (89.4 cm -1 ) for the B2 3 symmetry mode in the V point shows excellent 
agreement with the room temperature experimental data. This is to be expected since low energy 
modes show week anharmonicity effects. Surprisingly large discrepancy between experimental and 
calculated phonon energies is observed for the B3 5 mode. Since the B2 9 and B3 9 modes have similar 
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FIG. 4. (Color online) The Raman scattering spectra of FeSb2 (left panel) [15| and Feo.7sCoo.25Sb2 (right 
panel) single crystals in the (x'x>), x : = -^[101], configuration (A 9 symmetry modes) measured at various 
temperatures. 



normal modes, the chain rotation around the x and y axis, respectively (see Fig. [T]) their frequencies 
should be very close. This large disagreement is also unexpected since all other calculated phonon 
energies show rather good agreement with the experimentally obtained data. By detailed inspection 
of our previously published Raman spectra[l^] of pure, Co and Cr doped FeSb2 samples we did 
not find any mode in a low frequency region close to the calculated frequency (95 cm -1 ) for the 
mode. The missing B^ g mode is most probably of a very low intensity and it was not possible 
to extract it from the noise. The mode observed at 151.7 cm -1 for (x'y) polarization, which we 
assigned in Ref. Q as the B-^ g mode, could be the "leakage" of the A* mode, which appears at 
about 150.7 cm -1 in the (x'x 1 ) polarization. 

It is interesting to note that the dispersion curves of two A g symmetry Raman modes have 
opposites slopes near the T point (see the inset of Fig. E](b)), which leads to the mode mixing with 
the " anticrossing" effect. A* mode represents stretching vibration of Sb ions, whereas A^ mode 
represents twist ing of Sb ions which tent to rotate Sb ions around the z-axis, see Fig [TJ In our 



previous paper, 



15J we showed the existence of the A g mode mixing in the case of pure FeSb2 and 
Cr alloyed samples. Here we present detailed analysis of the mixing of two A g modes for pure and 
25% Co alloyed samples. 
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FIG. 5. (Color online) Energies (a), (b) and normalized intensities (c), (d) as a function of temperature of 
the A g modes for FeSb2 and Feo.7sCoo.25Sb2 single crystals. Thin solid lines show energy vs temperature 
dependence of the A g modes without coupling. Thick solid lines (red and black) show mode energy tem- 
perature dependance for two coupled A g modes calculated using Eq. The dashed lines are guide to the 
eye. 

The polarized Raman scattering spectra for pure FeSb2 (left panel [if]]) and Feo.7sCoo.25Sb2 (right 
panel) single crystals, measured in the (x'x') configuration (A g modes) at different temperatures, 
are presented in Fig. [H The Lorentzian lineshape profile has been used for the extraction of 
mode energy and linewidth. Fig. [5] shows the energies and normalized intensities as a function 
of temperature of the A g modes for FeSb2 and Feo.75Coo.25Sb2 single crystals. In the observed 
temperature range, energies of two A g modes for pure and 25% Co doped samples are very close 
which implies the existence of the mode mixing, manifested by mode repulsion and intensity transfer 
with the change of temperature. [28] Indeed, intensities of these modes are exchanged for both 
samples in the temperature range between 210 and 260 K (see Figs. 0]and[5]). 

In general, two phonon branches or any other elementary excitations of the same symmetry, may 
couple leading to the renormalization of the quasiparticle energies. Coupling between two phonon 
branches yields to the energy and linewidth changes (anticrossing effect). We can consider the 
coupling of two phonon branches as coupling of two quantum oscillators. When the perturbation 
is small, we can write the Hamiltonian of the system as 



H 



fii(T) V 

v n 2 (T) 



(3) 



where V is the interaction constant, fii(T) and ^(T) are the unperturbed mode energies, obtained 
by takeing into account, due to simplicity, only three-phonon process in Eq. (JT]). The eigenvalues 
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TABLE II. Best fit parameters for energy temperature dependance of the A g symmetry modes using Eq. 
©■ 



Compound 


Symmetry 


f2o (cm 1 ) 


C (cm" 1 ) 


V (cm- 1 ) 




A 1 


167.1 


2.65 




FeSb 2 








1.2 




A 2 


161.5 


1.16 






A 1 


167.5 


2.2 




Fe .75Co .25Sb 2 








2.4 




A 2 

A 9 


161.9 


0.80 





of the Hamiltonian are given by 

u± = i (fti (T) + Q 2 (T) ± V(^i(T)-^ 2 (T)) 2 + 4^ . ( 4 ) 

Eq. ([4]) gives a rather good fit of the experimental data (solid lines in Fig. E](a),(b)), suggesting 
the absence of any additional temperature dependent couplings (i.e. electron-phonon interaction) 
for these modes. Fitting parameters are presented in Table (TTJ Zero-temperature energies of A 1 
and Ag symmetry modes, in the absence of interaction, for pure (25% Co dopped) sample are 
fii(0) = 160.3 cm- 1 and n 2 (0) = 164.5 cm" 1 (Qi(0) = 161.1 cm" 1 and fi 2 (0) = 165.3 cm" 1 ). 
One can notice that the zero-temperature energies for decoupled modes are increased by 0.8 cm -1 
(about 0.5% increase) with 25% Co dopping, corresponding to the unit cell volume contraction. [3 1 
The phonon energy of the bond-stretching mode scales as R~ 3 , where R is the bond length. [29] 
Since the change in R~ 3 is proportional to the inverse volume change, we can expect the phonon- 
energy change for bond-stretching modes (A g modes) to be inversely proportional to the volume 
change. Because the Co atom substitutes Fe atom, which is located in the center of the inversion, 
there is no change in Raman spectra due to the mass effect. Additional repulsion between the 
coupled modes are due to the interaction. With Co doping, the interaction constant V increases, 
resulting in larger mode separation for the 25% dopped sample. 

V. CONCLUSION 

In summary, we presented a detailed theoretical and experimental study of the FeSb2 phonon 
dynamics. All experimentally observed Raman and infra-red active modes were successfully as- 
signed. The calculated phonon frequencies in the T point agree with the measured frequencies. We 
believe that the low energy B^g mode is of a very low intensity and therefore is not observed in the 
Raman experiments. The phonon mode at 150.7 cm -1 , which we previously assigned as the B^ g 
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mode, could be the "leakage" of the A* mode. The strong intensity exchange of the A g symmetry 
modes, observed in our Raman scattering experiments in the temperature range between 210 K 
and 260 K, is successfully described by a simple model of coupling of two phonon branches with 
the same symmetry. The mode mixing is also implied from the calculated dispersion curves, which 
show opposite slopes for two A g modes near the T point. We find that doping of FeSb2 with Co 
increases the A g modes repulsion. 
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